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Abstract 

A topological nodal-line semimetal is a new state of matter with one-dimensional bulk nodal lines and 
two-demensional “drumhead” surface bands. Based on first-principles calculations and our effective k ■ p 
model, we theoretically propose the existence of topological nodal-line fermions in the ternary transition- 
metal chalcogenide TlTaSe2. The noncentrosymmetric structure and strong spin-orbit coupling give rise to 
spinful nodal-line bulk states which are protected by a mirror reflection symmetry of this compound. This 
is remarkably distinguished from other proposed nodal-line semimetals such as Cu3NPb(Zn) in which the 
nodal line exist only in the limit of vanishing spin-orbit coupling and thus is not possible in real material. 
We show that the “drumhead” surface states in TlTaSe2, which are associated with the topological nodal 
lines, exhibit an unconventional chiral spin texture and an exotic Lifshitz transition as a consequence of 
the linkage among multiple “drumhead” surface-state pockets, suggesting a new type of 2 D topological 
metal. 
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Recently the experimental discoveries of three-dimensional (3D) topological Dirac semimetals 
and Weyl semimetals [l-20| have stimnlated enormous research interest in topological semimetals. 
Topological semimetallic (TS) materials are characterized by robust bulk band crossing points and 
the associated topological boundary states. In the 3D TS materials, the band crossing points can 
be either zero dimensional (OD) discrete nodal points or one dimensional (ID) continuous nodal 
lines. Materials that host these exotic band structures exhibit unique properties and hold promise 
for applications such as topological qubits, low-power electronics and spintronics. One prominent 
example of TS materials with OD band crossing points is the Weyl semimetals. The nodal points of 
Weyl semimetals carry non-zero chiral charges and are connected by the Fermi arc surface states. 
The Weyl semimetals have been experimentally realized in transition metal monopnictides such 


as TaAs 




14| . The other class of topological semimetals is the nodal-line semimetal whose 


conduction and valence bands cross each other at closed lines instead of discrete points. The 
topological nodal-line semimetal is distinct in three aspects compared to the Weyl semimatal: (1) 
the bulk fermi surface is ID and OD in nodal-line semimetals and Weyl semimetals, respectively; 
(2) the density of states of low-energy bulk excitations is proportional to {E — EfY and \E — Ef \ in 
nodal-line and Weyl semimetals, respectively; (3) The nodal lines are accompanied by “drumhead”- 


like surface states while Weyl nodal points are connected by ID Fermi arc surface states 2ll-l26|. 


The unique properties of nodal-line semimetals offer a new playground for studying novel physics 
arising from correlations between the massless quasiparticles. For example, interaction-induced 
instabilities that have been broadly discussed for Weyl semimetals should be more likely occurring 
in nodal-line states due to the higher density of states at the Fermi energy. 


Up to date, there have 


nodal-line semimetals 


Deen several theoretical proposals for a material realization of topological 


27 


31| . All these works predict nodal-line bulk states and “drumhead” 


surface states. However, the stability of nodal lines in the works requires the absence of spin-orbit 
coupling (SOC). With the inclusion of SOC, each nodal line is gapped due to the interaction 
between spin components. In real materials SOC, on the other hand, is ubiquitous, therefore it 
is important to study nodal-line semimetals under the condition of nonvanishing SOC. Generally, 


spinful nodal lines are unstable, which can be seen from a simple codimension analysis 


2^. In 


order to have robust nodal lines in the presence of SOC, an extra crystalline symmetry is needed 
to protect them. In this work, we report, based on hrst-principles calculations, the existence of 
spinful topological nodal lines in the ternary transition-metal chalcogenide TlTaSe 2 . The nodal 
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lines are 0.22 eV below the Fermi level and protected by a mirror reflection symmetry of the space 
group. The topological nodal-line state in TlTaSe 2 is in the class A+R (p = 2) of symmetry- 


protected semimetals 


23 


25| and is conneted with spin-polarized “drumhead” surface modes. 


We also demonstrate an exotic Lifshitz transition as a consequence of the linkage among multiple 
“drumhead” surface-state pockets. Our results establish TlTaSe 2 a promising material for studying 
nodal-line physics. 

TlTaSe 2 , shown in Fig. la, crystalizes in a hexagonal lattice in which the unit cell consists 
of one Tl, one Ta and two Se atoms and each atom resides on a hexagonal layer. The stacking 
sequence of these atomic planes within the unit cell is Tl-Se-Ta-Se: B-A-B-A (A, B and C, here, 
refer to the three high-symmetry spots on a hexagonal lattice). In other words, the Tl layer 
intercalates between adjacent TaSe 2 layers with Tl atoms aligned with Ta atoms in the vertical 
direction. The structure is non-centrosymmetric and belongs to the space group PQm2 (187). 
The lattice is reflection-symmetric with respect to both the Ta plane and the Tl plane. This 
reflection symmetry plays a key role in protecting the topological nodal lines, as discussed later 
on. Figure lb shows the bulk and (OOl)-projected surface Brillouin zones where the A, H and L 
points are high symmetry points on the plane, a mirror plane of the bulk Brillouin zone. 

The calculated band structure of TlTaSe 2 without and with the inclusion of SOC is shown in Fig. 
Ic. Around the H point, a hole pocket derived from Ta-5(ia,j^/a,2_y2 orbitals crosses an electron 
pocket from orbitals, taking the plane parallel to the Ta atomic plane as the x-y plane. 

All these atomic orbitals are invariant under the mirror reflection with respect to the Tl atomic 
plane. A zoom-in view of band structure around H is plotted in Fig. Id. In the case without 
SOC, the conduction and valence bands belong to different representations of the space group. A' 
and A" for electron and hole bands, respectively. The intersection of the two bands is, therefore, 
protected by the crystalline symmetry, forming a spinless nodal ring on the mirror plane kz = vr. 
Upon turning on SOC, each band splits into two spin-polarized branches since the system lacks 
space inversion symmetry. The spin splitting results in an accidental band touching of Ta and 
Tl bands at H and three band crossings. Only the band crossing 0.22 eV below the Fermi level 
remains gapless while the other two are gapped. A detailed analysis on the symmetry, orbital 
composition and spin texture of the bands around H is presented in Fig. 2a. The bands are mainly 
comprised of Ta.- 5 dxy/x 2 _y 2 and T\-Qpx^y states which are mostly confined the Ta and Tl atomic 
planes, respectively. The spin of these states is primarily oriented along z as indicated in Fig. 2a. 
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At the gapless crossing point, the two branches have opposite mirror parity eigenvalues. Therefore, 
this band crossing is protected by the mirror symmetry, forming a pair of nodal rings on the mirror 
plane kz = tt, one around H and the other around H'. To visualize the nodal rings, we plot the 
iso-energy contour aX E = —0.18 eV, shown in Fig. 2b. The energy is slightly off the nodal-line 
energy (—0.22 eV), which creates a tubular Fermi surface enclosing nodal lines. Indeed, we can 
hnd two rings surrounding H and H' points on kz = plane. At the center of the bulk Brillouin 
zone there is a spherical Fermi surface which is the hole pocket around F mainly derived from the 
Ta- 5 d 3 ;j 2_^2 orbitals. 

To further illustrate the mechanism of mirror-symmetry protection on the nodal rings, we 
develop an effective k ■ p description. The reflection symmetry Rz that sends to —in position 
space acts on the spin degree of freedom like the third Pauli matrix < 73 . The crystal structure of 
TlTaSe 2 is layered in the 2 : direction, where T1 and Ta layers are alternating. A consistent choice 
of the action of Rz is 


Rz'4^T\\nx,ny,nz ‘4^Tl\nx,ny,—nz-i Rz'4^Ta\nx,ny,nz '4’'Ta,\nx,ny,—nz — li (1) 

where ' 0 Ti|n^,ny,n 2 and i)Ta\n,:,ny,nz are the single-particle wave functions of an electrons in the T1 
and Ta orbitals in the unit cell indexed by the integers Ux, Uy, Uz- For Bloch states in the kz = n 
plane, the action of Rz is local in momentum space. If i^kz^XyX^ ^ 4-spinor in the space of spin 
and Tl/Ta orbital degrees of freedom, then Rz has the representation 

Rz'4^kx,kyj'^ TS ® (^) 


where the Pauli matrices cr* act on the spin space and the Pauli matrices Tj act on the Tl/Ta orbital 


space. 

We now give an effective k ■ p Hamiltonian to explain how the line nodes centered around the 
H/H' points arise. We will focus on H, implying that the Hamiltonian around H' follows as its time- 
reversal symmetric conjugate. We observe from the hrst principal calculations that the dispersions 
for the T1 and Ta orbitals are particle- and hole-like, respectively, around these high-symmetry 
points. In absence of spin-orbit coupling, the Hamiltonian for small momenta away from H is given 




kl + kl 


(To 4- vkzTi ® (To, 


which commutes with the representation 0 of Rz within the mirror plane. 
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We will include spin-orbit coupling terms that are constants to lowest order in k. Relevant 
mass terms are given by 

ky, k^) = Asoc'Ts ® a^ + Mi + M 2 , (4) 


where the mass terms Mi and M 2 can be any matrix that anticommutes with the kinetic term 
T 3 ( 8 ) (To, but commutes with the mirror symmetry ts 0 (T 3 . This leaves four choices of mass terms, 
namely Ti 0 (Ti, Ti 0 (T 2 , r 2 0 (Ti, and r 2 0 (T 2 . All of them will be in principle present. We have 
included two of them in the above Hamiltonian to reproduce the qualitative feature observed in 
the DFT calculations that the band gap above the nodal ring is much smaller than the band gap 
below the nodal ring. This means that two competing (i.e., commuting) mass terms are almost 
of the same magnitude. For example, we can choose Mi = miTi 0 ai and M 2 = 777 , 2 X 2 0 02 as 
competing mass terms. 

Let us consider for k^ = 0. The bands are given by 

I / k^ + k"^ ^ 

EniK, ky^ 0) = ±^soc ±Jimi + 7722)2 ^ - /ij , (5) 


with uncorrelated signs. The middle two bands are degenerate along (at most) two lines with 
kl + kl = 2m[fi ± v/A 2 q^ — (7771 -|- 777 , 2 ) 2 ]. These two nodal lines are protected by the opposite 
mirror eigenvalues of the crossing bands. The actual band structure of TlTaSe 2 is close to the 
parameter regime where the line node with smaller radius shrinks to zero. It disappears at the 
parameter point 


^soc — + ("^1 + "^ 2 )^- (6) 

The band structure from the effective Hamiltonian is plotted in Fig. 2c, which is in remarkable 
agreement with the hrst-principles result. Considering the crystalline symmetry of TlTaSe 2 , the 
nodal line at H is classihed as the time-reversal breaking class A+R which admits an integer 
topological classihcation for Fermi surfaces of codimension 2, i.e., lines (p = 2 in 23|, l2^). We 
note that although the entire system preserves time-reversal symmetry, the effective Hamiltonian 
at H lacks time-reversal symmetry, which makes the nodal lines of TlTaSe 2 in a different class 
from those reported in 28|, |30(|. The nodal lines of TlTaSe 2 are characterized by a topological 
quantum number 77 ^, which is given by the difference in the number of occupied bands with 
eigenvalue -|-1 inside and outside the nodal line. In the case at hand, 77 ^ = -|-1 for the nodal line 
at H in the k^ = 'n' plane as shown in Fig. 2c. We also perform a band structure calculation in 
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which that the Ta atom is slightly moved in the vertical direction and, thus, the mirror reflection 
symmetry is broken. In this case all of the four branches around H are found to belong to the 
same S2 representation of the reduced space group and a gap opening is allowed at every crossing 
point of these branches, as shown in Fig. 2d. Therefore, the nodal line in this case is gapped by 
mirror-symmetry-breaking perturbations. In other words, the nodal rings in TlTaSe 2 are indeed 
under the protection of the mirror reflection symmetry. 

Next we study the evolution of the nodal lines as SOC varies. The nodal-line band structure 
at various SOC is plotted in Fig. 3a. As shown in the above discussion, without SOC there is a 
spinless nodal line around H. Once turning on SOC, each band becomes two spin branches with 
opposite mirror parity eigenvalues and the spinless crossing line is, consequently, split into four 
spinful crossing lines. However, among these four crossing lines only two are robust under the 
protection of mirror reflection symmetry. In this case they are the crossings between bands bl 
and b4 and between b2 and b3, because the mirror parity of b2 and b4 is -|-1 and that of bl and 
b3 is —1. So there are two nodal rings around H arising from the two symmetry-protected band 
crossings. The crossings between bl and b3 and between b2 and b4 are gapped. The gap size at 
the topmost bl-b3 crossing is very small, but it can be evidently seen in Fig. 2a and explained by 
our effective Hamiltonian as shown in Fig. 2c. As SOC increases, bands b2 and b3 are gradually 
pulled apart from each other and with SOC = 1 in the scale relative to the real SOC of the material 
the two bands barely touch. For even larger SOC, e.g. SOC = 1.2, the two bands separate and, as 
a result, the nodal ring associated with these two bands disappears. Therefore only one nodal ring 
is left when SOC is beyond the critical value 1. In Fig. 3b, we plot the energy difference between 
the top of band b2 and the bottom of band b3 as a function of SOC. Varying SOC, we have three 
different phases of nodal lines, namely, spinless nodal line, double nodal line, and single nodal line. 

In order to illustrate connection of “drumhead” surface states with the bulk nodal lines, the sur¬ 
face electronic structure is constructed using a flrst-principles-derived tight-binding model Hamil¬ 
ton in a slab. The projected (OOl)-bulk band structure along T — K — M direction is plotted in 
Fig. 4a where the nodal line is clearly seen. The surface can be terminated with either a Se layer 
or a T1 layer, because the bonding between Se and T1 atoms is much weaker than that between 
Se and Ta atoms. The surface band structure corresponding to the two possible terminations is 
shown in Figs. 4b and 4c. In the case of Se termination, the surface band disperses outwards with 
respect to K from the nodal line and grazes outwards at the edge of the bulk band below the nodal 
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line and merges into the bulk band region, as marked by the arrows in Fig. 4b. This can be also 
seen in the iso-energy band contour at E = —0.25 eV (slightly below the energy of the nodal line, 
-0.22 eV), shown in Fig. 4d. The “drumhead” surface states overlap with the bulk nodal rings and 
there is no other surface state at this energy. On the other hand, the surface band structure on 
the Tl-terminated surface is dramatically different. Along T — K direction, starting from the nodal 
point, the surface band grazes outwards at the edge of the upper bulk Dirac cone and merges into 
the bulk band. Along K — M direction, the “drumhead” surface band SSI disperses into the bulk 
band gap. Within the gap there exists a second surface band SS2 which joins the “drumhead” 
band at M forming a Kramers pair. The spin polarization of the two surface bands is shown in Fig. 


4e. The two surface bands possess opposite spin polarizations. The overall surface banc 


and spin texture resemble those of the Dirac surface states of topological insulators 


32 


structure 


33j. The 


“drumhead” surface band is without spin degenracy and it arises from the band inversion of two 
spinful bulk bands (bl and b4). Figure 4f shows the iso-energy band contours of Tl-terminated 
surface at three different energies (indicated by dotted lines in Fig. 4e). At E = —0.25 eV, slightly 
below the energy of the nodal rings, two surface bands form a human-eye-shaped contour, lying in 
between two nodal-ring bulk pockets at K and K'. The spin texture is unique in the sense that as 
moving along either surface band clockwisely, the spin orientation rotate in counterclockwise di¬ 
rection. The spin orientation is not always in the tangential direction of the Fermi surface contour, 
which is distinct from the spin-momentum-locked texture of Dirac surface states in conventional 
topological insulators |32|, l33|. At the energy of the nodal rings {E = —0.22 eV), the eye-shaped 
contour connects to the nodal rings and, as energy moves further up, the corners of the eyes open 
and the “drumhead” surface band SSI transforms into a giant closed contour surrounding f. This 
topological change in the band contour is known as Lifshitz transition in the electronic structure. 
The Lifshitz transition discussed here is special in two aspects: (1) the transition and the associate 
saddle-point singularity happen in spin-polarized 2D surface bands; (2) the transition relies on a 
linkage of multiple small pockets around M to form a giant pocket surrounding F. This exotic link¬ 
age of 2D surface bands give rise to a divergence in DOS, like the case in the van Hove singularity in 


topological crystalline insulators and high-Tc superconductors 


34 


37| . Shifting chemical potential 


to the Lifshitz transition energy in TlTaSe 2 by means of chemical doping or electrical gating can 
potentially trigger interaction-induced instabilities such as unconventional superconductivity 361 
or spin/charge density waves at the surface. 
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In summary, topological nodal-line semimetals is a distinct class of topological materials beyond 
topological insulators and Weyl semimetals. In this work we propose theoretically that TlTaSe 2 , a 
ternary transition-metal chalcogenide, is a promising candidate for material realization of topologi¬ 
cal nodal-line semimetals. Unlike previous proposed materials, the nodal lines in TlTaSe 2 is robust 
even with the inclusion of SOC as long as the mirror reflection symmetry with respect to the Ta 
atomic plane is not broken. The nodal line is 0.22 eV below the Fermi level and, thus, accessible 
to the conventional ARPES measurements. Through systematic surface simulations, we show the 
unique spinful “drumhead” surface states on Se-terminated surface and, even more interestingly, 
hnd that Tl-terminated surface of TlTaSe 2 features a Lifshitz transition as a consequence of the 
exotic linkage of multiple “drumhead” surface-state pockets. In light of these novel properties 
of the electronic band structure of TlTaSe 2 , we establish an ideal material platform for studying 
unique nodal physics and electronic correlation effect of topological nodal-line semimetals. 


METHODS 


We computed electronic structures using the norm-conserving pseudopotentials as implemented 


in the OpenMX package within the generalized gradient approximation (GGA) schemes 38 


Experimental lattice constants were used 


401 1 . Al2xl2x4 Monkhorst-Pack k-point mesh was 


used in the computations. The SOG effects are included self-consistently [^. For each T1 atom, 
three, three, three, and two optimized radial functions were allocated for the s, p, d, and / orbitals 
{s3p3d3f2), respectively, with a cutoff radius of 8 Bohr. For each Ta atom, d3p2d2fl was adopted 
with a cutoff radius of 7 Bohr. For each Se atom, d3p2d2fl was adopted with a cutoff radius of 7 
Bohr. A regular mesh of 300 Ry in real space was used for the numerical integrations and for the 
solution of the Poisson equation. To calculate the surface electronic structures, we constructed hrst- 
principles tight-binding model Hamilton. The tight-binding model matrix elements are calculated 


by projecting onto the Wannier orbitals 


42|. We use T1 p, Ta s and d, and Se p orbitals were 


constructed without performing the procedure for maximizing localization. 


ACKNOWLEDGEMENTS 

T.R.G. acknowledges visiting scientist support from Princeton University. We thank Ghuang- 
Han Hsu for technical assistance in the theoretical calculations. We thank Ghen Fang and Andreas 








10 


P. Schnyder for discussions. 


[1] Weyl, H. Elektron und gravitation I. Z. Phys. 56, 330-352 (1929). 

[2] Volovik, G. E. The Universe in a Helium Droplet (Oxford University Press, 2003). 

[3] Herring, C. Accidental Degeneracy in the Energy Bands of Crystals. Phys. Rev. 52 , 365 (1937). 

[4] Wilczek, F. Why are there analogies between condensed matter and particle theory? Phys. Today 
51 , 11-13 (1998). 

[5] Balents, L. Weyl electrons kiss. Physics 4, 36 (2011). 

[6] Hasan, M. Z., Xu, S.-Y. &: Neupane, M. in Topological Insulators: Fundamentals and Perspectives 
(eds Ortmann, F., Roche, S. &: Valenzuela, S. O.) (John Wiley &: Sons, 2015). 

[7] Wan, X., Turner, A. M., Vishwanath, A. &: Savrasov, S. Y. Topological Semimetal and Fermi-arc 
surface states in the electronic structure of pyrochlore iridates. Phys. Rev. B 83 , 205101 (2011). 

[8] Burkov, A. A. &: Balents, L. Weyl semimetal in a topological insulator multilayer. Phys. Rev. Lett. 
107,127205 (2011). 

[9] Xu, G. et al, Chern semi-metal and quantized anomalous Hall effect in HgCr 2 Se 4 . Phys. Rev. Lett. 
107, 186806 (2011). 

[10] Singh, B., Sharma, A., Lin, H., Hasan, M. Z., Prasad, R. Sz Bansil, A. Topological electronic structure 
and Weyl semimetal in the TlBiSe 2 class of semiconductors. Phys. Rev. B 86 , 115208 (2012). 

[11] Huang, S.-M., Xu, S.-Y. et al, An inversion breaking Weyl semimetal state in the TaAs material 
class. Nature Commun. 6:7373 (2015). 

[12] Xu, S.-Y. et al., Discovery of a Weyl Fermion semimetal and topological Fermi arcs. Science 349, 
613 (2015). 

[13] Weng, H., Fang, C., Fang, Z., Bernevig, A. &: Dai, X. Weyl semimetal phase in non-centrosymmetric 
transition metal monophosphides. Phys. Rev. X 5 , 011029 (2015). 

[14] Lv, B. Q. et al., Experimental Discovery of Weyl Semimetal TaAs. Phys. Rev. X 5 , 031013 (2015). 

[15] Xu, S.-Y. et al. Observation of Fermi arc surface states in a topological metal. Science 347, 294-298 
(2015). 

[16] Yang, B.-J. & Nagaosa, N. Glassification of stable three-dimensional Dirac semimetals with nontrivial 
topology. Nat. Commun. 5, 4898 (2014). 



11 


[17] Wang, Z., Weng, H., Wu, Q., Dai, X. &: Fang, Z. Three-dimensional Dirac semimetal and quantum 
transport in Cd 3 As 2 . Phys. Rev. B 88, 125427 (2013). 

[18] Neupane, M.et al., Observation of a three-dimensional topological Dirac semimetal phase in high- 
mobility Cd 3 As 2 . Nat. Commun. 5, 3786 (2014). 

[19] Borisenko, S. et al., Experimental realization of a three-dimensional Dirac semimetal. Phys. Rev. 
Lett. 113 , 027603 (2014). 

[20] Liu, Z. K. et al., Discovery of a three-dimensional topological Dirac semimetal, Na3Bi. Science 343, 
864-867 (2014). 

[21] Burkov, A. A., Hook, M. D. & Balents, L. Topological nodal semimetals. Phys. Rev. B 84 235126 

( 2011 ). 

[22] Phillips M. & Aji V. Tunable line node semimetals. Phys. Rev. B 90, 115111 (2014). 

[23] Chiu, C. &: Schnyder, A., Classification of reflection symmetry protected topological semimetals and 
nodal superconductors. Phy. Rev. B 90, 205136 (2014). 

[24] Fang, C., Chen, Y., Kee, H.-Y. &: Fu. L. Topological nodal line semimetals with and without spin- 
orbital coupling Chen. arXiv:1506.03449v2. 

[25] Bian, G. et al.. Topological nodal-line fermions in the non-centrosymmetric superconductor PbTaSe2. 
arXiv: 1505.03069 (2015). 

[26] Yang, S. A., Pan, H., & Zhang, F. Dirac and Weyl Superconductors in Three Dimensions. Phys. Rev. 
Lett. 113 , 046401 (2014). 

[27] Xie, L. S., Schoop, L. M., Seibel, E. M., Gibson, Q. D., Xie, W. &: Cava R. J. Potential ring of Dirac 
nodes in a new polymorph of Ca 3 P 2 . APL Materials 3 , 083602 (2015). 

[28] Yu, R., Weng, H., Fang, Z., Dai X. & Hu X. Topological Node-Line Semimetal and Dirac Semimetal 
State in Antiperovskite Cu 3 PdN. Phys. Rev. Lett. 115 , 036807 (2015). 

[29] Zeng, M., Fang, C., Chang, G., Ghen, Y.-A., Hsieh, T., Bansil, A., Lin, H. & Fu, L. Topological 
semimetals and topological insulators in rare earth monopnictides. arXiv: 1504.03492 (2015). 

[30] Kim, Y., Wieder, B. J., Kane, C. L. & Rappe, A. M. Dirac line nodes in inversion-symmetric crystals. 
Phys. Rev. Lett. 115 , 036806 (2015). 

[31] Weng, H. et al.,Topological node-line semimetal in three-dimensional graphene networks. Phys. Rev. 
B 92, 045108 (2015). 

[32] Hasan, M. Z. &: Kane, C. L. Topological insulators. Rev. Mod. Phys. 82, 3045-3067 (2010). 


12 


[33] Qi, X.-L. & Zhang, S.-C. Topological insulators and superconductors. Rev. Mod. Phys. 83, 1057-1110 

( 2011 ). 

[34] Okada, Y. et al., Observation of Dirac node formation and mass acquisition in a topological crystalline 
insulator. Science 341, 1496-1499 (2013). 

[35] Hirsch. J.E., Scalapino. D. J., Enhanced Superconductivity in Quasi Two-Dimensional Systems. 
Phys. Rev. Lett. 56, 2732-2735 (1986). 

[36] King. D. M. et al, Observation of a Saddle-Point Singularity in Bi2(Sro.97Pro.03)2Cu06+5 and Its 
Implications for Normal and Superconducting State Properties. Phys. Rev. Lett. 73, 3298 (1994). 

[37] Liu. C. et al, Evidence for a Lifshitz transition in electron-doped iron arsenic superconductors at 
the onset of superconductivity. Nature Phys. 6, 419 (2010). 

[38] Perdew, J. P. Burke, K. &: Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. 
Rev. Lett. 77, 3865 (1996). 

[39] Ozaki, T. Variationally optimized atomic orbitals for largescale electronic structures. Phys. Rev. B 
67, 155108 (2003). 

[40] Eppinga, R. &: Wiegers, G. A. A generalized schene for niobium and tatalum dichalcogenides inter¬ 
calated with pos-transition elements. Physica B and C 99, 121 (1980). 

[41] Theurich, G. Sz Hill, N. A. Self-consistent treatment of spin-orbit coupling in solids using relativistic 
fully separable ab initio pseudopotentials. Phys. Rev. B 64, 073106 (2001). 

[42] Weng, H., Ozaki, T. & Terakura, K. Revisiting magnetic coupling in transition-metal-benzene com¬ 
plexes with maximally localized Wannier functions. Phys. Rev. B 79, 235118 (2009). 



13 



FIG. 1: Crystal lattice and bulk band structure of TlTaSe 2 . (a) Lattice structure of TlTaSe 2 : top 
view (left) and side view (right), (b) Bulk and (OOl)-projected surface Brillouin zones, (c) Calculated 
bulk band structure of TlTaSe 2 without (left) and with (right) the inclusion of spin-orbit coupling. The 
color shows the atomic orbital decomposition, (d) Zoom-in band structure around H point as marked in 


c. 
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FIG. 2: Symmetry, orbital composition and spin textnre of nodal-line Fermions in TlTaSe 2 . 

(a) Orbital composition, space-group representation and spin polarization of nodal-line bands in TlTaSe 2 . 
The mirror parity of each band is given in parentheses, (b) Iso-energy bulk-band contour at E = -0.18 
eV of TlTaSe 2 bulk bands, (c) Band structure of the effective k ■ p model given by Eqs. m and f|4|) that 
approximates the low energy bands around the H point for parameters mi = 0.7, m 2 = 0.8, /i = 2.2, 
^SOC = 2.5, with Rz reflection eigenvalues of the bands and reflection-protected crossings (orange circles) 
that corresponds to the nodal line. All bands are nondegenerate and the momentum k = -|- ky, kz = 0 

is measured relative to the high-symmetry point H. (d) Calculated bulk band structure of TlTaSe 2 with 
the Ta atom shifted slightly away from the equilibrium position in the unit cell. The shift breaks the 
mirror symmetry of the system and reduces the symmetry of the nodal-line states. 
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AE=E(b3)-E(b2) 



FIG. 3: Evolution of the bulk band structure of TlTaSe 2 as the spin-orbit coupling increases. 

(a) The nodal-line band structure of TlTaSe 2 with various SOC. (b) The energy difference between the 
two bands b2 and b3 as a function of SOC, showing nodal-line phases at different SOC. 
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FIG. 4: Nodal-line fermions and “drnmhead” snrface states of TlTaSe 2 . (a) (OOl)-projected bulk 
bands of TlTaSe 2 . (b) Bulk and surface band structure of Se-terminated surface of TlTaSe 2 . The surface 
bands are marked by the arrows, (c) Same as b bnt for Tl-terminated snrface of TlTaSe 2 . (d) Iso-energy 
band contour at E = -0.25 eV of Se-terminated surface. The “drumhead” surface states overlap with the 
bulk nodal rings, (e) The spin polarization of surface bands of Tl-terminated surface, (f) Iso-energy band 
contour at E = -0.25 eV (top), -0.22 eV (middle), and -0.18 eV (bottom) of Tl-terminated surface. The 
arrows indicate the spin texture of the surface bands. 

















